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Oh' Abstract. In this article the NA61/SHINE detector will be presented which is a large 

acceptance hadron spectrometer at the CERN SPS. It allows for a precise study of particle 
production from interactions of a 30 GeV proton beam in a carbon target in order to predict the 
neutrino flux of the T2K experiment at J-PARC, Japan. Requirements for the T2K experiment 
will be discussed together with the ongoing NA61 measurements. In particular preliminary 
^ ' NA61 results on pion production and cross section measurements from the 2007 pilot run will 

C , be shown. 
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The NA61/SHINE (SHINE = SPS Heavy Ion and Neutrino Experiment) experiment at the 
CERN SPS combines a rich physics program in various fields [TJ [2j [31 [3] . Besides performing 
measurements for the T2K experiment, it takes a variety of data used for the description of 
cosmic-ray air showers in the Pierre Auger and KASCADE experiments [5j [6] as well as for 
studying the behaviour of strongly interacting matter at high density. This article focuses 
on the NA61 hadron production measurements in p+C interactions which are targeted for an 
accurate neutrino flux prediction of the T2K experiment [7]. 

T2K is a long baseline neutrino oscillation experiment at J-PARC, Japan, with the aim to 
precisely measure the v„ — > u e appearance and disappearance. To generate neutrinos a high 
intensity 30 GeV proton beam impinging on a 90 cm long carbon target is used, whereby mesons 
(it, K) are produced which decay into neutrinos {v^ e )- The neutrino flux is then measured 2.5° 
off-axis both in a near detector, located 280 m behind the target, and the 295 km far away Super- 
Kamiokande detector (SK). Neutrino oscillations can be probed by comparing the neutrino flux 
measured at SK to the one predicted at SK. In order to predict the flux at SK one uses the near 
detector measurements and extrapolates them with the help of Monte Carlo (MC) predictions 
to SK. Up to now, these MC predictions depend on hadron production models. For more precise 
predictions measurements of pion and kaon production off the carbon target are essential. The 
aim of the NA61/SHINE measurements for T2K is to provide this information. 

The NA61 detector is a large acceptance hadron spectrometer at the CERN SPS (see 
Fig- !])• It's main components originate from the NA49 experiment [SJ. These include four 
large volume time projection chambers (TPCs), which provide high momentum resolution 
( a (p)/p 2 — (0-3 — 7) • 10~ 4 (GeV/c) -1 ) and good particle identification (PID) via specific 
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Figure 1. The layout of the NA61/SHINE set-up (top view, not in scale). 

energy loss measurements. The PID capability of the TPCs is complemented through velocity 
measurements of two time-of-flight (T0F-L/-R) detector arrays with a time measurement 
resolution of a ~ 60 ps. In 2007 an additional ToF detector (ToF-F) with a resolution of 
a ~ 120 ps was installed in order to extend the PID acceptance to low momenta (p < 4 GeV/c) 
as required for T2K |4|]. Two carbon (isotropic graphite, p = 1.84 g/cm 3 ) targets of different 
lengths are used: 



• a 2 cm long target (~4% of nuclear interaction length, A/), the so called thin target, and 

• a 90 cm long cylinder of 2.6 cm diameter (~1.9 A/), the so called T2K replica target. 

Proton beam particles are identified and selected by means of CERDAR and threshold Cherenkov 
counters as well as several scintillation counters. Interactions in the thin target are selected by 
an anti-coindidence of the incoming beam particle with a small scintillation counter (S4) of 2 cm 
diameter placed on the beam axis between the two vertex magnets. 

During the pilot run [3] in October 2007 approximately 660k events with the thin target 
(target in), 80k events without target (target out) and 220k events with the T2K replica target 
were registered with a 30 GeV proton beam. From the observed interaction probabilities for 
target in (Pr in ) and out (-Pr otlt ) configurations a trigger cross section (atrig) has been evaluated 
as: 

°tri 9 = pLef l fNA / A " P Tou t ) = (297.5 ± 0.7(stat.)± < 10%(syst.$) mb , (1) 

where Na, p, A and L e ff denote, respectively, the Avogadro constant, the target density, it's 
atomic number and it's effective length. The inelastic cross section (ai ne i) can be obtained 
from a^ig by applying two major corrections accounting for trigger biases. The first one is 
to subtract the elastic contribution due to large angle coherent scattering of primary protons 
which do not reach the trigger counter S4 even though no inelastic interaction occured. The 
second one is to add the lost inelastic contribution coming from secondary particles hitting S4 
and therefore preventing from triggering on the event. Both corrections have been estimated 
with GEANT4 simulations and result in Oi ne i = (254.7 ± 1.0(stat.)± < 10% (syst.) 1 ) mb. atrig 
and ai ne i are finally used for the normalization of the particle spectra. Up to now, preliminary 
spectra for positively and negatively charged pions have been obtained via three different particle 
identification methods based on: 



1 A more precise evaluation of the systematic error is under investigation. 
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Figure 2. Upper panel: dE/dx as a function of momentum for p+C at 31 GeV/c beam 
momentum for positively and negatively charged particles, respectively, together with the Bethe- 
Bloch curves. Lower panel: Examples of combined ToF and dE/dx PID performance for 
positively charged particles in the momentum range 2-3 GeV/c (left), 3-4 GeV/c (middle) and 
4-5 GeV/c (right). 

• Energy loss measurements in the TPC gas (see Fig. [2]). A dedicated dE/dx analysis for 
particles with momentum below 1 GeV/c, which is the region not covered by the time-of- 
fhght detectors, was performed. 

• Combined energy loss and time-of- flight measurement information (dE/dx + ToF). 
It was used to perform identification of pions with momenta above 1 GeV/c (see Fig. [2]). 

• Analysis of negatively charged particles further referred to as h~~ . It is based on the 
theoretical and experimental premises that negative particles produced by 30 GeV protons 
consist mainly of negative pion mesons with a small admixture of electrons, negative kaons 
and a negligible fraction of antiprotons. This procedure allows to obtain spectra of it~ 
mesons in the full momentum range. 

For the three different PID methods a common Venus-GHEISHA and Geant Monte Carlo 
simulation chain is used to calculate corrections for geometric acceptance, reconstruction 
efficiency, weak decays and vertex association. Preliminary results on tt + and ir~ differential 
inclusive inelastic cross sections obtained from the different analyses are presented in Figs. [3] 
and HI respectively. Only statistical errors are shown. The spectra obtained using the different 
methods agree better than by 20%, which yields the upper limit for systematic biases. Work 
to minimize these biases is currently in progress. Measurements obtained with the T2K replica 
target are currently being analysed. In 2009 further 6 million interaction triggers for the thin 
target and 4 million for the T2K replica target were recorded, which will increase the limited 
statistics from the 2007 run by a factor of about 10. 
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Figure 3. Differential inclusive inelastic cross sections for tt + versus momentum in different 
angle intervals from dE/dx (pink) and dE/dx + ToF (blue) analyses. 
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Figure 4. Differential inclusive inelastic cross sections for it versus momentum in different 8 
angle intervals from dE/dx (red) and dE/dx + ToF (green) and h~ (black) analyses. 
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